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Abstract. In this paper, we simulate the linear absorption spectra of the 
MePTCDI and PTCDA crystals. The basic Hamiltonian describes the Frenkel 
excitons and charge-transfer excitons mixing in the molecular stack (point group 
Ci) and their linear coupling with one vibrational mode of an intramolecular vi- 
bration. Using the vibronic approach, we calculate the linear optical susceptibil- 
ity in the excitonic and one-phonon vibronic regions of the molecular stack and 
of a crystal with two types of nonequivalent stacks. We put the excitonic and 
vibrational parameters for the two crystals fitted in the previous studies and ana- 
lyze some important features of of the linear absorption lineshape in the spectral 
region of 15 000-23 000 cm^^ and the virtual positions of many-particle bands. 
Our study exhibits the necessity of introducing the FE-CTEs mixing in the in- 
terpretation of the linear absorption spectra, especially in the MePTCDI crystal. 

PACS codes: 71 .35.Aa, 73.20.Mf, 78.40.Me 



1 Introduction 

The mixing of Frenkel excitons (FEs) and charge-transfer excitons (CTEs) has 
been studied both experimentally and theoretically in many one-component 
molecular crystals, e.g., polyacenes |[l}21, perylene derivatives ||5j^2l, fullerenes 
|l8]|9l, and other. The excitonic and vibronic spectra of quasi-one-dimensional 
crystals like 3,4,9, 10-perylenetetracarboxylic dianhydride (PTCDA) and N, N'- 
dimethylperylene-3,4,9,10-perylenetetracarboximide (MePTCDI) have been 
treated in and |fT0HT4l . The short molecular distance in quasi-one- 

dimensional stacks causes a strong FE-CTEs mixing and mixing of their vi- 
bronic satellites as well, especially in the absorption spectra in the spectral re- 
gion of 2-3 eV (15 000-23 000 cm^^). 



*This work is dedicated to the memory of Academician Matey Mateev (1940-2010), a great 
friend and wonderful colleague on Bulgarian physics community, 
t Corresponding author: izh@phys . uni-sof ia . bg 
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In the present paper, we calculate the linear absorption spectra of one-component 
MePTCDI and PTCDA crystals applying the vibronic approach developed in 
our previous papers II13II151 . It allows complex calculations of the absorption 
in pure excitonic, one-phonon vibronic, two-phonon vibronic, etc. spectra. In 
paper 1161 the vibronic approach is the tool of studying the FE-CTEs mixing in a 
two-component stack of alternatively arranged donor-acceptor (DA)-molecules. 
While the FE-CTEs mixing in fTSl is a probable but hypothetical model, in the 
present study we turn to the real absorption spectra and use the parameters of 
the FE-CTEs mixing from Refs. lISllelfTOl. 

Two general differences exhibit the models of the FE-CTEs mixing in one- and 
two-component molecular stacks ITSl . notably: (i) The intermolecular transfer 
both of electrons and holes must be a feature of the model for the case of one- 
component stack whereas only one transfer mechanism is sufficient in a DA- 
molecular stack. In the last case the Frenkel exciton represents a collectivized 
electronic excitation of a donor or an acceptor and the strongest transfer on the 
closest neighbor would be either of a hole or correspondingly of an electron. 
Obviously in an one-component stack both types of transfer on the neighbors 
could be of the comparable probability. The two-step processes of successive 
transfer of the electron and the hole on the same molecule ensure the transfer of 
a EE even in the case of its relatively weak direct intermolecular transfer ifTOl . In 
this way, our present study extends the calculations in ifTJl where only one trans- 
fer mechanism has been considered, (ii) In the most widely studied DA-crystal, 
anthracene-PMDA the excitonic absorption lines are narrow lfT7tjT9l and many 
details in the vibronic spectra can be seen. In PTCDA and MePTCDI crystals 
the absorption lines are two order of magnitude wider Thus we pay attention to 
the general structure of the excitonic and one-phonon vibronic spectra suppos- 
ing a width of the excitonic lines of 300-500 cm^^ (not 2-10 cm^^ which is the 
absorption width of the anthracene-PMDA). 

In the next section of the paper we involve the initial Hamiltonian in the case 
of a FE-CTEs mixing. The Hamiltonian contains one mode of intramolecular 
vibration linearly coupled to the FE and CTEs. In Section 3 the linear opti- 
cal susceptibility has been calculated in the excitonic and one-phonon vibronic 
spectra. In Section 4 the linear absorption has been modelled using the exci- 
tonic and vibrational parameters of the PTCDA and MePTCDI crystals fitted in 
Refs. 161171 [TOl . Section 5 summarizes the main findings derived in this paper 
In the Appendix we treat some problems of the FE-CTEs mixing in a molecular 
stack of the Ci point group of symmetry. 

2 Hamiltonian for thie Case of a FE-CTEs-Phionon Coupling 

We consider the excitonic and vibronic excitations in a linear molecular stack 
of identical molecules which are regularly arranged at a distance d each 
other The point group of symmetry of the stack is Ci as is for the molecular 
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stacks of PTCDA and MePTCDI f^. The origin of the Frenkel exciton is a 
non-degenerate molecular electronic excitation with excitation energy E-p and 
L being the transfer integral between neighboring molecules. We denote by £?„ 
(5+) the annihilation (creation) operator of the electronic excitation on molecule 
n and get the following FE-part of the Hamiltonian: 

HfE = ^ EyB^Bn + ^ L ((5„'.„+l + <5n',„-l) B^,Bn. (1) 
n n.n' 

As usually, we consider two CTEs of equal excitation energy Ec and C„.i is the 
annihilation operator CTE, a = 1, with hole located on the site n and electron 
on the site n + 1, whereas the electron of the second CTE, cr = 2, (C„.2) is 
located at molecule n — 1. We neglect the transfer of CTEs as a whole and the 
mutual coupling of the two CTEs since those processes can be realized through 
the transfer of the electron or hole at distance 2d which is less probable than the 
FE-CTEs mixing caused by the electron (hole) transfer at the neighbor molecule 
(see Refs. II16II18I ). We have the following CTEs-part of the Hamiltonian: 

-ffcTE = ^ EcC:^crCna (2) 
n, 17=1,2 

and suppose the following operator for the FE-CTEs mixing: 

n 

+ e^iB+Cn-1,1 + eh2S+C„+i,2 + h.c] , (3) 

where Sci and ec2 are the transfer integrals of the electron from molecule n to 
molecules n + 1 and n — 1 correspondingly, and Shi, £h2 denote the transfer 
integrals of the hole from molecule n to molecules ?? + 1 and n — 1. Certainly 
the model with four transfer integrals is more complicated than the model in 
Refs. miTl with two mixing parameters only. But our model is more realistic 
because we take into account the inclination of the flat molecules of PTCDA and 
MePTCDI relatively to the stack axis (see the Appendix). 

One intramolecular mode is only supposed to be coupled with the EE and CTEs 
and the phonon part of the Hamiltonian is 

Hph = ^ /iwoa+a„, (4) 

n 

where ujq is the vibrational frequency and a„ is the annihilation operator of one 
vibrational quantum on molecule n. The linear exciton-phonon coupling only is 
manifested in the treated crystals Q and we get the following exciton-phonon 
part lUiTlEol 

^ex-phon = ^ fit^o [(.pB^Bn (a+ + a„) 

n,a=l,2 

+ CC'naCna (a+ + a„ + a+_^^^ + ttn+ai)] , (5) 
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where <^f and ^ are dimensionless parameters of the hnear FE-phonon and 
CTEs-phonon coupUng, correspondingly; cti = +1 if cr = 1, and cti = — 1 
if (7 = 2. For the sake of simplicity, we suppose the same linear exciton-phonon 
coupling in the positive and negative ions which create CTEs. 

The full Hamiltonian H contains all the parts ([B-© and can be transformed 
using the canonical transformation which eliminates the linear exciton-phonon 
coupHng Q, see II20II21I . 

Hi=cxp{Q)Hcxp{~Q), (6) 

where 

Q= [epB+Sn (a+ - a„) 

+ ^C+^Cna (a+ - an + 0++,^^ - a„+^J] . (7) 

We introduce the vibronic operators 

K =exp(g)B„exp(-g), (8) 

Un^a = exp{Q)C„cr exp{-Q) (9) 
and get the following transformed Hamiltonian 

Hl^J2 ~ ^"^f) VuVn {Sr,,n+l + (5„-,„_i) V+Vn 

n nn' 

n,fT n n 

+ £c2V;+t/„2 + ehiK^ t/„-i,i + £h2K^ t/„+i,2 + h.c] (10) 

In a stack with inversion center, point group Ci, the excitons are gerade or unger- 
ade (in the center of the Brillouin zone, at fc = 0). The ungerade excitons only 
are dipole-active and influence the linear optical susceptibility and the absorp- 
tion spectra. In the case under consideration, the operator of the transition dipole 
moment has the following form II16II17I : 

^ = E + Vn) + PCT (f/„2 - C/„l + (7+ - (7+ )] . (11) 

n 

The gerade EEs, as well as the symmetrical combination of CTEs, (U^2 
+ U^-^ |0), can be also mixed, but due to their vanishing transition dipole mo- 
ment they will not be considered here. 

Introducing the Eourier transform in the momentum space of the vibronic oper- 
ators, 

yk = ^Y. ^" cxp(ifcnd), (12) 
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Uk.a = -7= ^ Un,a exp(ifcnd), (13) 



we obtain the following form of the Hamiltonian (fTOl l 

Hi=Y^ {Ef - ftwoCi + cos k) V+Vk 

k 

k.a 

+ 5' + ^ricjoa+a„, (14) 

n 

where 

S^Y. {^^^^ (^^-1 + + BV+ {Uk,2 - Uk,i) + h.c.} 

with 

A — e'f, + e'l^cosk + ie'^sink and -B = Eg + e(( cos fc + iej^ sin /c. 

Here 

4 = (Sci +£c2)/2, £h = (eiii +eh2)/2, (15) 
£c = (£o2 - £cl) /2, £h = (£h2 - £hi) /2. (16) 

For the case of ungerade FEs, their mixing with the symmetrical (even) combi- 
nation {Uk,i + Uk;2) at fc = is impossible, and thus the mixing parameters 
and vanish, = = 0. The final expressions for operators (|3} and (ITll are 

^FCT = [«' + COSfc) V+ {Uk,2 - Uk,l) 

k 

+ ie'^ sin kV+ {Uk,2 + Uk,i) + h.c] (17) 

and 

P = ViV [pF (Ffe=o + 14to) 
+ PCT (c/fc=o,2 - C/fe=o,l + C/+^o,2 - C^fcto,i)] - (18) 

respectively. 

3 Calculation of the Linear Optical Susceptibility 

The linear optical susceptibility can be calculated by using the formula ll22l 

= l™ I ^ + ie) + $y + ie)] I (19) 

5 



/. J. Lalov, I. Zhelyazkov 



with 

^,j{t) = -ie{t){o\p,{t)Pj{o) + Pj{t)P,{o)\o), (20) 

where V is the crystal's volume [in our case being proportional to Nv {v is 
the volume occupied by one molecule)] and P is the operator ([18). The Green 
functions ( |20] | have been calculated as average over the ground state |0) only by 
taking into account the large values of Ep, Ec, hoJo ^ kBT. 

We calculate the Green functions (l20l i following the vibronic approach II15II16I . 
In the next expression, the x axis is supposed to be oriented along the vector pp 
which includes angle 7 with the vector pcT, and a = |pct/pf|- Then we can 
represent the linear optical susceptibility of one stack as 

Xxx = - — ["2 + 4aai2 COS 7 + 2a^ai cos^ 7] . (21) 

The functions ai, a2, 0112 have been calculated for the excitonic and one-phonon 
vibronic regions (see below). The PTCDA and MePTCDI crystals contain two 
types A and B of parallel molecular stacks, however, the excitonic and vibronic 
excitations in each stack interact very weakly with the excitations of the other 
stacks. In the same way as in Ref. lfT3l . we calculate the crystal's susceptibihty 
in an oriented gas model. We denote by 2ip the angle between the vectors pp 
and pf of two different stacks and suppose that these vectors are positioned in 
the [XY) plane, the crystal X axis been oriented along the sum pp + pf \22>\. 

The components of the linear optical susceptibility of the crystal correspond- 
ingly are: 

XXX = :t-y- "2 cos if 

V aia2 — 2ai2 

+ 4aai2 cos^ ip cos 7 + a^ai (1 + cos 27 cos (p)] (22) 

and 

Pf 2 r • 2 

XYY = [a2 sm ip 

+ 4aai2 sin^ (y9C0S7 + a^ai(l — cos27COS(p)] . (23) 

We find the following expressions for functions ai, a2, ayi- 

(1) In the excitonic region expressions practically coincide with formulas in 
Ref. HSl, namely 

ai =fia;-(i;F + 2L)-nr2oF(l), (24) 

«i2=£e+eh, (25) 
a2 - ft[w-l^oc(l)] -i^c, (26) 
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where f2oF(l) and r2oc(l) are expressed through the continuous fractions fol- 
lowing from recursions: 



Lo - (£:f + 2L) /h - nujo - ^oFin + 1) ' 
^^0c(?^) = 



Lu — Ec/h — nujf) — ^oc{n + 1) ' 

Wa = ^F^O, Wal = ^Wq. 



(2) In the one-phonon vibronic region 



ai = huj -{Ey + 2L) 



mD 



[hujic - I3)a 



I \ 'Mil ~a"-ai I V 



(27) 

(28) 
(29) 

(30) 
(31) 



a2 = /zw - i?c - ^ (— - QF'T + [e'l + elf 
mD [cjic 

cr + (Ti(2 + <) a-pCFi 



where 



= w - Wo - r2ic(i) - Ec/h, 



(32) 



(33) 



m = 2 [Lhwic + 2e^'ej;] 



(34) 



(1/to) \{h [lj-lo^- r!iF(l)] - -Bf} hioi, - 2 [(4')' + Kf] } , (35) 



D= 1- 



1 - 



2 r 



apCi (7 + (7i(2 + t) 



(36) 



and 



-i/\/F~T if k + Vt2~~T|<i, 



i/\/F~T if k + \/<2~~T| < 1, 



(Tl = to" — 1. 



(37) 



(38) 
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Table 1. Excitonic and vibrational parameters (in cm ^) of the MePTCDI and PTCDA 
crystals. 





Ep 






L 




Eh 


2lp 


7 


MePTCDI 
PTCDA 


17 992 

18 860 


17 346 

18 300 


1400 
1400 


345 
330 


-380 
-48 


-137 
-436 


36.8° 
82° 


68.1° 
143° 



The functions r^iF(l) and riic(l) also represent continuous fractions from re- 
cursions: 

ullJ 

uj — hp /n — (n + ljujQ — i^lF(?^ + Ij 

u! — hc/n — [n + l)ujo — ihc{n + Ij 
Finally, the functions /3 and ap can be expressed as follows; 

/3/ft=f}oc(2)-r!ic(l), (41) 
aFA=r!oF(2)-r!iF(l). (42) 

4 Simulations of the Excitonic and Vibronic Spectra of lUlePTCDI 
and PTCDA Crystals 

In this section, we calculate the absorption spectra of the two crystals finding 
the imaginary parts of the components of the linear optical susceptibility ( |22] | at 
(2pp/v) = 1 and supposing an imaginary part equal to i6/h of the frequency uj. 
We put the excitonic and vibrational parameters for the studied crystals as they 
have been fitted in Refs. ^ and ifTOI and used in our previous papers lfTTHT4l . 
see table 1 : The data for angle 7 and ratio a = \pcT /pf \ have been calculated 
in 161 using quantum chemical evaluations. We use the values a ~ 0.1 for the 
PTCDA crystal and a = 0.135 for the MePTCDI. The data for angle 2ip are 
derived from the crystal structure (for MePTCDI see ll24l '). 

The aforementioned parameters are permanent in our calculations. We vary the 
values of the following parameters (intending to observe their impact on the 
absorption spectra and find a better similarity with the experimental absorption 
spectra lEJCl): 

(i) The excitonic damping quantity S which varies from 1 to 500 cm^^. 

(ii) The linear exciton-phonon coupling parameters and ^ which vary near 
the values 

= 0.82 and 1.1 for PTCDA (see 12), 
= 0.88 and 1.1 for MePTCDI (see Q), 



8 



Excitons-CTEs-Phonons Coupling in . 



and 

^ = ^F, S,f/V2, and ^f/2 for both crystals. 
In calculating continuous fractions ( |27| i. ( [28] l and (|39] l. (|40t we take twenty steps. 

4.1 MePTCDI 

The fine structure of the excitonic spectra of this crystal is presented in Figure 1 . 
The absorption curves are calculated at small value 6 = 20 cm~^ of the exci- 
tonic damping for = 0.88 with ^ = (blue curve 1), ^ = ^f/\/2 (red 
curve 2) and ^ = (green curve 3). We have calculated the imaginary part 
of expression (l22l i using the 'excitonic' formulas (I24l -(|26]| (in the following de- 
noted as "exc" program). For the smaller values of ^ equal to / \/2 and ^f/2 
we obtain spectral doublets which correspond to the CTEs-FE splitting. This 
splitting generates a spectral triplet at 15 100 cm^^ (very weak), 16 500 cm^^, 
and 17 200 cm^^ (for ^ = ^p, blue curve 1). 

0.025 
0.02 

0.01 
0.005 


15000 15500 16000 16500 17000 17500 

CO (cm"^) 

Figure 1. Linear absorption in the excitonic region of the MePTCDI crystal at 5 = 20 
cm"\ = 0.88 with C = Cf (blue curve 1), ^ = ^f/V^ (red curve 2) and ^ = 
(green curve 3) calculated using formulas (24)-(26) ("exc" program). 

The value 5 = 300 cm^^ seems to be more closer to the width of the lowest 
absorption maximum observed in the MePTCDI crystal ||6][2]. The absorption 
curves in Figure 2 are calculated at 5 — 300 cm^^ in the excitonic region (ap- 
proximately below 17800 cm~^), one-phonon vibronic region (17 800-19 200 
cm~^), two-phonon vibronic region (19 200-20 700 cm~^) etc. There exists an 
obvious similarity between experimental curve (see Ref. ^) and the green curve 
3 in Figure 2 calculated with £, ~ /2. 

The impact of the CTEs transition dipole moment is illustrated in Figure 3. The 
two curves - the red one 1 calculated with the CTEs contribution, and the green 
one 2 calculated without this contribution [a = in formulas (l2ni-(l23]ll - are 
relatively close to each other The strongest impact can be measured near lower 
excitonic peak generated by the CTEs level = 17 346 cm^^. 
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0.0025 
0.002 
^ 0.0015 
0.001 
0.0005 


15000 16000 17000 18000 19000 20000 21000 22000 23000 

CO (cm"^) 

Figure 2. Linear absorption in the excitonic region of the MePTCDI crystal at 5 = 300 
cm~^ ("exc" program). For the values of and ^ see Figure[T] 

0.0025 
0.002 
^ 0.0015 
0.001 
0.0005 


15000 16000 17000 18000 19000 20000 21000 22000 23000 

CO (cm"^) 

Figure 3. Impact of the CTEs-transition dipole moment expressed through the ratio a. 
Red curve 1 for a = 0.135 and the green curve 2 for a = (5 = 300 cm^^, — 0.88, 

e = ^f/2). 

The two curves in Figure 4 calculated by using the "exc" program for S = 300 
cm~^ differ by the mutual positions of the two excitonic levels. The green curve 
1 corresponds to a lower position of the CTEs level Ec < Ep, whereas the 
red curve 2 corresponds to the opposite situation E-p < Ec (the magnitudes of 
Ec and Ep are exchanged). The red curve 2 exhibits a strong domination of the 
lower maxima associated with the FE level. The two curves can be approximated 
with five Lorentz maxima in the spectral region of 15 000-22 000 cm~^ 16]. 
Comparing the lineshapes of these two curves with the experimental curve in 
Ref. 161, we cannot make a hypothesis which possibility is more probable. We 
prefer the fitting from Refs. ||6] and IfTOl {Ec < Ep), however, another choice, 
Ep < Ec, is also allowed. 

Figure 5 illustrates the impact of the FE-phonon coupling constant on the 
absorption spectra. The values = 0.88 and ^ — £,p/2 can be considered as 
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0.0025 
0.002 
^ 0.0015 
0.001 
0.0005 


15000 16000 17000 18000 19000 20000 21000 22000 23000 

CO (cm"^) 

Figure 4. Comparison of the absorption curves calculated with the data from Table 1 
(green curve 1) and the red curve 2 calculated for = 17 346 cm^^ and Ec = 17992 
cm~^ (interchange of the excitonic levels) with 5 — 300 cm~^, ~ 0.88, ^ = ^f/2. 

0.002 

0.0015 

^ 0.001 
E 

0.0005 


15000 16000 17000 18000 19000 20000 21000 22000 23000 

CO (cm"^) 

Figure 5. The impact of the FE-phonon coupling constant on the linear absorption. 
The red curve 1 coiTesponds to ^f = 0.88, and the green curve 2 to ^f = 1-1 (^ = Cf/2, 
5 = 300 cm"^). 

better candidates for the simulation (compared with the experimental absorption 
lineshape especially in the excitonic region). 

The importance of the FE-CTEs mixing can be seen in Figure 6 in which the 
blue curve 2 represents a pure FE-absorption (e" = e[j = 0). The complexity 
of the experimental absorption curves with five Lorentz maxima in the studied 
spectral region, see Figures 2 and 3 in Ref. [i6j|, cannot be understood on the basis 
of simple Frenkel-exciton model. Contrary, the model of mixed FE-CTEs and 
their vibronic satellites can be the basis of adequate simulations of the absorption 
curves (red curve 1). 

The absorption curves in the following three Figures 7-9 have been calculated 
by using formulas (|27]|-(|42]| for one-phonon vibronic spectra denoted in the fol- 
lowing as "Ip" program. Our calculations based on the Green functions formal- 
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0.003 



0.0015 - 

E 

0.001 - 



0.0005 



0.0025 - 



0.002 



" 




15000 16000 17000 18000 19000 20000 21000 22000 23000 




Figure 6. Linear absorption spectra of the MePTCDI crystal with FE-CTEs mixing (red 
curve 1) and pure Frenkel-exciton spectra (blue curve 2 at e" ~ £h ~ 0) obtained by the 
"exc" program with 5 = 300 cm~\ = 0.88, C = ?f/2. 



Figure 7. One-phonon vibronic spectra ("Ip" program) of the MePTCDI crystal with 
— 0.88, ^ = ^f/\/2, S — 1 cm~^ (red curve 1), 5 = 2 cm~^ (green curve 2), and 
5 = 300 cm~^ (purple curve 3). 

ism allow to study and simulate the absorption associated with the two types of 
exciton-phonon states in the vibronic spectra, notably: (a) bound states (one- 
particle states) corresponding to the propagation of excitons and phonon(s) in 
the stack as a whole. The lineshape of their absorption maxima is Lorentzian 
and strongly depends on the excitonic damping 6. (b) Many-particle states (MP) 
corresponding to the excitation of exciton and phonon(s) on separate molecules 
of the stack. Their energy lies in the quasicontinuous band(s) and the lineshape 
of the corresponding absorption maxima depends on the density of the states in 
the MP band and more weakly on the excitonic damping 6. Unfortunately the 
big values of 6 in both crystals under consideration mask the MP bands which 
are totally covered by the absorption maxima. We can simulate the absorption 
in MP continua by supposing very small values of the damping, for instance. 




0.0025 







17400 17600 17800 18000 18200 18400 18600 18800 

CO (cm"^) 
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S ^ 1-10 cm ^, which is not typical for the perylene derivatives (but occurring 
in some DA-crystals Hke An-PMDA imi2l[T7l[T8ll). 



0.007 
0.006 
0.005 
0.004 
0.003 
0.002 
0.001 



A1 

M 


(a) ■ 





19400 19500 19600 19700 19800 19900 20000 20100 20200 



CO (cm"^) 




CO (cm"^) 



Figure 8. (a) MP band in two-phonon vibronic spectra with — 0.88, ^ = ^f/\/2, 
(5 = 1 cm~^ (red curve 1), 5 = 10 cm^^ (blue curve 2). (b) Bound exciton-phonon state 
in one-phonon vibronic spectra with ~ 0.88, — £,f/ V2, 3 = 1 cm^^ (red curve 1), 
3 = 2 cm^^ (green curve 2), 5 = 10 cm^^ (blue curve 3). 




18000 19000 20000 21000 22000 23000 

CO (cm"^) 

Figure 9. Linear absorption in the vibronic spectra of the MePTCDI crystal ("Ip" pro- 
gram) with FE-CTEs mixing (red curve 1) and pure FE vibronics (green curve 2 at 
e" = = 0) with S = 300 cm-\ Cf = 0.88, ^ = ?f/2. 



The one-phonon vibronic spectra are depicted in Figure 7. The vibronics of FE 
and CTEs are mixed hke their excitons' spectra. Nevertheless, the less intensive 
lower maximum at 17 500-17 700 cm~^ can be considered as vibronic replica 
of CTEs whereas the more intensive maximum near 18 400-18 700 cm^^ is the 
replica of the FE. The lower wings of both maxima represent the very weak MP 
bands (even virtual ones). The main part of absorption intensity is concentrated 
above those MP bands in the Lorentzian maxima near 17 600 cm~^ and 18 650 
cm^^. As can be seen in Figure 7, the wide absorption maxima corresponding 
to the realistic value S = 300 cm^^ (purple curve 3) cover the fine structure of 
the absorption lineshape. 
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Figure 8 illustrates the difference between the MP band (left panel (a)) and the 
bound exciton-phonon maximum (right panel (b)). The line shape of the MP 
band depends of the density of the states which is strongly modulated by the 
exciton-phonon coupling. Thus, the lineshape is non-symmetrical. In the case of 
relatively small damping parameter S the lineshape (look at Figure 8(a)) depends 
very weakly on 6. Contrary, the maximal value of the Lorentzian maximum in 
Figure 8(b) is proportional to 1/6. 

Figure 9 is analog of Figure 6 in the vibronic region (above 17 500 cm~^). The 
lineshape of the red curve 1 (with FE-CTEs mixing) is more complicated than 
the lineshape of the pure FE vibronics. The red curve 1 seems to be a better 
simulation of the experimental absorption curve ||6]. 

The model of the FE-CTEs mixing reproduces the general structure of the ab- 
sorption spectra in the excitonic and vibronic regions of the MePTCDI crys- 
tal (see Figures 2 and 9). Our calculations confirm the correct choice of the 
excitonic and vibrational parameters in the fitting procedure implemented in 
Refs. 161 [TOl . especially the value = 0.88 of the constant of the linear 
FE-phonon coupling. Additionally, we establish as the most probable value 
of the exciton damping 6 = 300 cm^^ and CTEs-phonon coupling constant 
^ = Our studies show the positions of the bands both of FE+phonon and 

CTEs+phonon many-particle states. However, the wide absorption lines cover 
the virtual MP bands. 

4.2 PTCDA 

The experimental linear absorption lines for the PTCDA crystal are broader than 
those of the MePTCDI O. Thus, the structure of the PTCDA spectra does not 
exhibit many details generated by the FE-CTEs mixing. 

Our calculations show the appearance of the spectral doublets in the excitonic 
and vibronic spectra associated with the FE-CTEs splitting (red curve 1 in Fig- 
ure 10). However, the probable value of the excitonic damping parameter is 
6 = 500 cm^^ evaluated using the width of the lowest experimental absorption 
maximum ||6l. The blue curve 2 in Figure 10 calculated with 6 = 500 cm~^ ex- 
hibits four Lorentz-type maxima as it is observed in the experimental absorption 
spectra of the PTCDA crystal |6j. The fifth maximum near 17000 cm^^ (see 
Figure 10) manifests itself as a weak pecuUarity only of the absorption curve 
with S = 500 cm"^ 

Moreover, further calculations implemented with = 1-1 (which value seems 
to be the most suitable to our simulations) demonstrate insensitivity of the ab- 
sorption spectra on the value of ^ (see the very close green curve 1 and blue 
curve 2 in Figure 11). We also calculated the absorption curve neglecting the 
FE-CTEs mixing (purple curve 3). The result consists of some re-distributions 
of the absorption intensity among four maxima but the absorption lineshape is 
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Figure 10. Linear absorption of the PTCDA crystal calculated using the "exc" program 
with = 0.82, 5 = ^f/2, 5 = 10 cm"^ (red curve 1), and S = 500 cm"^ (blue curve 
2). 
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Figure 11. Linear absorption of the PTCDA crystal ("exc" program) with = 1-1, 
S = 500 cm^^, ^ = i^f/2 (green curve 1) and ^ = ^f (blue curve 2). The purple curve 
3 corresponds to pure FE spectra (e" — £]" — 0). 

practically non-changed in the cases with and without FE-CTEs mixing. 

The same effect - the lineshape being not strongly affected by a hypothetical 
FE-CTEs mixing - can be seen in Figure 12 (green curve 1 and blue curve 2 
there). The mutual replacement of the two excitonic levels (compare the red 
and green curves) shifts the position of the absorption maxima, however, their 
Uneshapes are very similar 

The absorption curves in Figures 13 and 14 have been calculated using the "Ip" 
program and they show the positions and the structure of one-phonon and two- 
phonon vibronic spectra. Surely, the fine structure can be observed in the hy- 
pothetical case of small values of S. Then the structure of the MP bands de- 
pends on the FE-CTEs mixing. Figure 13(a) containing two absorption regions, 
18 600-19 200 cm-i and 19 400-19 900 cm~^, each one with a weak MP band 
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Figure 12. Linear absorption of the PTCDA crystal ("exc" program) with — 1.1, 
^ = ^f/2, and 5 — 500 cm^^. The green curve 1 corresponds to Ef = 17 992 cm^^, 
Ec = 17 346 cm~^, the blue curve 2 to pure FE spectra (e" = ej' — 0), and the red 
curve 3 to = 17346 cm"\ Ec = 17992 cm~\ 
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Figure 13. One-phonon vibronic spectra of the PTCDA crystal ("Ip" program) with (a) 
^P = 0.82, ^ = ^p/2, S — 1 cm^^ (red curve 1), S = 10 cm~^ (green curve 2) and 
(b) = 0.82, ^ = ^f/2, and = e'^ — (vibronic spectra of FE). The red curve 1 
corresponds to (5 = 1 cm~^, the blue curve 2 to <5 = 10 cm^^, and the purple curve 3 to 

S = 500 cm"\ 

and a Lorentzian. This picture is rather different in the case of the absence of 
FE-CTEs mixing (Figure 13(b)). Then we might observe a CTEs vibronic max- 
imum near 18 800 cm^^ and the vibronic replica of FE with MP band (18 900- 
19 400 cm^^), as well as a Lorentz-type maximum at 19 600 cm^^. In our 
model Hamiltonian (|2]i the dispersion of pure CTEs has been neglected and thus 
at e" = e[[ = the vibronic replica of CTEs consists of one-particle maximum 
only. In the same time, EEs possess dispersion expressed through the transfer 
terms with L in expression ([U (see also (ISITOl) and in the vibronic absorption 
spectra of FEs both one-particle maxima and MP band could appear. As in other 
cases, a strong damping with 6 = 500 cm~^ masks the whole structure. 

The absorption curves in Figure 14 correspond to a higher value of = 1.1. 
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Figure 14. One-phonon and two-phonon vibronic spectra of the PTCDA crystal ("Ip" 
program) with = 1.1, ^ — ^f/2. The blue curve 1 corresponds to 5 = 1 cm~^ and 
the purple curve 2 to <5 = 500 cm^^. 

The one-phonon vibronic spectra (18 500-19 000 cm^^) and two-phonon vi- 
bronic spectra (19 700-20 500 cm^i) contain MP bands (18 300-18 500 cm^i 
and 19 700-20 000 cm~^) and several other maxima. Their splitting is caused 
by the relatively big value of the constant of the linear FE-phonon coupUng. 
In the realistic case of strong damping, 5 ~ 500 cm^^, the one-phonon and 
two-phonon vibronic spectra would contain two broad and asymmetric maxima 
only. 

The main conclusion of our simulations of the hnear absorption in the PTCDA 
crystal is that the hypothesis of a FE-CTEs mixing is not crucial in describing 
the observed spectra. The wide absorption lines in the excitonic and vibronic 
spectra cover the effects of the mixing and the hypothesis of Frenkel excitons' 
spectra and the vibronic of FE could be sufficient (see Ref. Il25l ). 

5 Conclusion 

Our model for the linear absorption spectra of the one-component charge- 
transfer molecular crystals includes the following parameters: 

• excitonic levels E-p and Ec, as well as the vibrational frequency Wq of the 
intramolecular mode, 

• parameters e" and ej^' of the FE-CTEs mixing, 

• constants and ^ of the linear exciton-phonon coupling, 

• the width 5 of the excitonic linewidth, 

• angles 2Lp, 7 and the ratio a of the CTEs and FE transition dipole mo- 
ments. 

Practically all parameters have been introduced and fitted for MePTCDI and 
PTCDA crystals in previous papers ||5l-|7][T0l. based primarily on the matrix di- 
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agonalization method. In our study we apply the complex vibronic approach 
based on the Green functions formalism in calculating the linear optical sus- 
ceptibility and its imaginary part which is a factor in the absorption coefficient 
in the excitonic and one-phonon vibronic spectra. Higher vibronics - with 
two, three phonons - also have been demonstrated in our calculations. Cal- 
culated vibronic spectra consist of one-particle Lorentzian maxima and many- 
particle (MP) bands which correspond to unbound propagation of the excitons 
and phonons. Our simulations expose the positions of the MP continua but the 
wide excitonic and vibronic absorption lines in both crystals cover the fine struc- 
ture of the vibronic spectra. The main goal of our model is to simulate the line- 
shape in the absorption region of 15 000-23 000 cm^^ of the aforementioned 
crystals, and to find out better fitted values of the excitonic linewidth, exciton- 
phonon coupling parameters and so on. 

We stress again the main results which concern the two crystals: 

1. The excitonic linewidth of the MePTCDI crystal, according to our simu- 
lations, is approximately S w 300 cm^^. The FE-phonon linear coupling 
coefficient has been estimated correctly in the previous papers li5j|7] 
as ~ 0.88 and ^ = Cf/2. The absorption spectra depend strongly on 
the mutual position of the two excitonic levels Ep and Ec- A supposition 
for the FE-CTEs mixing is necessary for an adequate interpretation of the 
excitonic and vibronic spectra of the MePTCDI crystal. 

2. The excitonic linewidth of the PTCDA crystal can be evaluated as (5 « 500 
cm^^. That is why the effect of the FE-CTEs mixing, being covered 
by the wide absorption maxima, are more weakly expressed than in the 
MePTCDI crystal. However, a very probable conclusion from our calcu- 
lations may be the stronger linear exciton-phonon coupling (^p ~ 1 or 
1.1 instead of 0.82). 

Our model can be applied in the interpretation of other one-component molec- 
ular stacks (crystals). It can be more effective in the systems with narrower ex- 
citonic absorption lines where both types of vibronic states - one-particle states 
and MP continua - calculated by using the Green functions formalism would be 
seen in the linear absorption. 

Appendix 

In this Appendix, we treat two groups of problems which concern the mixing of 
Frenkel and charge-transfer excitons in an one-component stack of the Ci point 
group of symmetry. 

(A) The first group is rather formal. In the Ci point group, at the center of 
Brillouin zone, fc = 0, electronic excitations of the stack possess the following 
characteristics; 
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1 . The symmetry of inversion center makes the parity of the excitations odd 
or even and produces ungerade and gerade excitons. In another part of 
the Brillouin zone these two types of excitons can mix producing some 
interesting effects (see, for example, Refs. ||261 and ifTSI ). 

2. The ungerade excitons only are dipole active and possess non-vanishing 
transition dipole moment. For CTEs the ungerade wave function is the 
combination imiUfTTl: 

3. The mixing of two excitons of different types of FE and CTEs is allowed 
only in the case of the same parity of the coupling excitons. 

Thus, the ungerade EEs as well as the ungerade CTEs only can mix and 
can be active in the linear absorption. In the same time, the gerade FE and 
the gerade combination of CTEs 

4 [Kk=, + t^l^=o) |0> = 4= E + Kn) |0). (A.3) 

can mix but they will be non-active in the linear absorption. 

(B) The second group of problems concerns the mixing constants Ed, £02. ^hi. 
and eii2 in a molecular stack of the Ci point group (see Eqs. (|3]l and (fTOl i in the 
text). 

Let us consider three neighbour molecules A, B, and D located at the positions 
n — l,n, and n + 1 of the stack: An-i, Bn, Dn+i- We denote by i?* a neutral 
excited molecule (with FE on it) and by A^, B^, - the wave functions of the 
corresponding positive and negative ions. The matrix elements of the electron 
transfer can be represented as follows: 

{An-iB:Dn+i\H\An-iB+D-+,) 

= {An-lB:^Dn+l\H\An-lCn,l) = Sol, (A.4) 

{A„^^B:D„+,\H\A-_,B+D^+^) 
= (A„_iS,*A.+i|i?|C„,2A.+i> = £c2, (A.5) 

where H is the Hamiltonian, C„ ^ are the wave functions of the CTEs, £oi and 
£o2 are the mixing constants (similar expressions hold for £hi and £h2)- 

The gerade FE (£?* is an even function relatively to the symmetry of inversion) 
mixes with the combination (l/\/2) (C,ti + ^^2) |0) ^'^'^ '^he mixing constant 
is 

^(£el+£e2). (A.6) 
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In the case of ungerade FE the mixing constant with the ungerade combination 

{1/V2) (C+2 - C+i) |0) is: 

^ (£c2 - eel) . (A.7) 

The manifestation of the FE-CTEs mixing depends on the conditions at which 
the following equalities hold: 

£oi = £o2 and £hi = £vi- (A.8) 

The disposition of the molecules in the stack allows two possible configurations 
in which relations (IA.8b are fulfilled: 

1) The planes of the molecules are perpendicular to the stack's axis. In that 
case the configurations B^D'^j^^ and A^_^B^^ are mirror-symmetrical 
with respect to the positive ion i3+ and hence the quantities ( IA.4I ) and 
(IA.5b are equal, independently on the symmetry of the molecular orbitals 
in the positive ion and, on the other hand, in the negative ion. Obviously 
the perpendicular position of the molecules in the stack is not the case of 
the treated crystals [25]. 

2) In the case of inclination of the molecules' planes relatively to the stack's 
axis the full symmetry of the positive ion _B+, governed by the symmetry 
of the molecular orbitals on it, must coincide with the full symmetry of 
the negative ions and -D,7+i (governed by their molecular orbitals). 
Especially the maxima of the charge distribution in the positive and nega- 
tive ions must be located at the same points of the entity. In the opposite 
case, because of the inclination, the conditions of the electron/hole trans- 
fer from the neutral excited molecule n on the its two neighbours n — \ 
and 71 + 1 are non-equivalent. 

Since situation 2) seems to be less probable than the non-coinciding symmetry 
of the positive and negative ions, the inequalities £ei 7^ £e2 and £hi 7^ £h2 hold. 
This is notably the case of mixing FEs and CTEs which are active in the linear 
absorption, too. In that case, formula (flTi describes the FE-CTEs mixing of the 
dipole active excitons. 

In the same time, the following operator describes the FE-CTEs mixing of ger- 
ade excitons (in all the configurations): 

^FC = [(^o + Eh COS fc) 1/+ (C/fc,i + C/fc,2) 

k 

+ i£;, sin fcy+ (C/fc,2 - Uk,i) + h.c] (A.9) 
(see formula (flSl l for the expressions of e'^ and e'^. 
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The spectra of the gerade mixing FE-CTEs are similar to the spectra calculated 
in this paper, but they will not manifest themselves in the linear absorption. 

If configurations 1) or 2) are realized, then relations (IA.8t are fulfilled and op- 
erator (IA.9) describes the FE-CTEs mixing among gerade excitons only. In 
this case, the spectra of the linear absorption consist of two independent (non- 
mixing) spectra, notably: the linear absorption spectra of FEs and their vibron- 
ics, and the spectra of the CTEs and their vibronics. The simulation of non- 
mixing FEs and CTEs can be performed by using the formulas in the present 
paper, in which we have to put e" — e'^ = 0. 
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